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We consider decay of metastable states of forced vibrations of a quantum oscillator close to 
bifurcation points, where dissipation becomes effectively strong. We show that decay occurs via 
quantum activation over an effective barrier. The decay probability W scales with the distance 
to the bifurcation point as | In W| oc rj^ . The exponent ^ is found for a resonantly driven oscillator 
and an oscillator modulated at nearly twice its eigenfrequency. 
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Decay of a metastable state is usually considered as 
resulting from tunneling or thermal activation. In this 
paper we study a different decay mechanism, quantum 
activation. It relates to systems far from thermal equi- 
librium. As tunneling, quantum activation is due to 
quantum fluctuations, but as thermal activation, it in- 
volves diffusion over an effective barrier separating the 
metastable state. 

Metastable decay in noncquilibrium systems has at- 
tracted much attention recently in the context of switch- 
ing between coexisting states of forced vibrations. Such 
diverse systems as trapped electrons and atoms P, 0, 
Josephson junctions 0,^, and nano- and micromechani- 
cal oscillators 0, have been studied. The experiments 
largely focused on the parameter range where the system 
was close to a bifurcation point in which the metastable 
state disappears. In this range the decay probability is 
comparatively large and displays characteristic scaling 
with the distance to the bifurcation point. So far classi- 
cal activation was studied, but recently quantum regime 
has been also reached 0- 

For classical systems, scaling of the rate of activated 
decay near a bifurcation point was found theoretically 
both in the cases of equilibrium |^, 0, 0| and nonequi- 
librium systems ll3 • the latter case a scaling 

crossover may occur as the system goes from the under- 
damped to overdamped regime while approaching the bi- 
furcation point |l4j . Such crossover occurs also for quan- 
tum tunneling in equilibrium dissipative systems [l5j . 

In this paper we study decay of metastable vibrational 
states in dissipative systems close to bifurcation points, 
where the motion becomes overdamped. The analysis 
refers to the systems of current interest, quantum oscil- 
lators driven by a resonant force or parametrically mod- 
ulated at nearly twice the eigenfrequency. We show that 
at low temperatures decay occurs via quantum activa- 
tion. The decay rate W scales with the distance to the 
bifurcation point rj as \ lnW\ (x rj^ . The scaling exponent 
is ^ = 3/2 for resonant driving, and ^ = 2 for parametric 
modulation; in addition, jlnM^j displays a characteristic 
temperature dependence. 

Quantum activation in periodically modulated systems 
can be understood by noting that metastable states are 



formed as a result of the balance between external driv- 
ing and dissipation due to coupling to a thermal bath. 
For T = dissipation corresponds to transitions to lower 
energy states with emission of excitations of the bath. 
However, modulated systems are more adequately de- 
scribed by the Floquet (quasienergy) states than by the 
energy eigenstates. Emission of bath excitations may re- 
sult in transitions to both higher and lower quasiener- 
gies, albeit with different probabilities The 
higher-probability transitions lead to relaxation towards 
a metastable state, whereas the lower-probability tran- 
sitions lead to effective diffusion away from it, a finite- 
width distribution over quasienergy, and metastable de- 
cay. There is certain similarity here with the Unruh ef- 
fect 01 where a uniformly accelerated relativistic detec- 
tor coupled to a quantum zero-temperature field is de- 
scribed in its proper time by the Gibbs distribution with 
the acceleration-dependent temperature. 

We will start with a resonantly driven nonlinear oscil- 
lator. Its Hamiltonian is 
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In the presence of weak damping the oscillator may have 
two coexisting stable states of classical forced vibrations 
jl9l | . They emerge already for a small modulation ampli- 
tude A provided the detuning dcu = cup — loq of the mod- 
ulation frequency LOp from the oscillator eigenfrequency 
luq is small, \Suj\ <C iop. We assume that the nonlinear- 
ity is small, |7|((?^) ^ ojq, and that -ySuj > 0, which is 
necessary for the onset of bistability. 

It is convenient to switch from q,p to slowly varying op- 
erators Q, P, using a transformation q = Crcs(Q cosujpt+ 
Psmujpt), p = —Crcs^p{Q sinajpt — P cosujpt) with 
CrcB = {8ujpSu;/3-/y^^. The variables Q, P are the scaled 
coordinate and momentum in the rotating frame. 



[P, Q] = -iX, A = ah-y/Sujp 6uj. 
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The parameter A plays the role of the effective Planck 
constant. We are interested in the semiclassical case; A 
is the small parameter of the theory, A <C 1. 

In the rotating wave approximation the Hamiltonian 
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for Slli > becomes Hq = {h/X)5uj g, with 

9^9{Q,P) = \{Q' + P'-lf-f3'/'Q, (3) 

(for (Jti) < one should redefine g —g, Hq 
— {h/X)6u!g). The function g plays the role of the oscil- 
lator Hamiltonian in dimensionless time r — t\Suj\. The 
eigenvalues of g give oscillator quasienergies. 

The parameter /? in Eq. is the scaled intensity of the 
driving field. For weak damping the oscillator is bistable 
provided < (3 < 4/27. In this range the function 
g{Q, P) has a shape of a tilted Mexican hat. The max- 
imum at the top of the central dome and the minimum 
at the lowest point of the rim correspond, respectively, 
to the small- and large-amplitude states of forced vibra- 
tions. The saddle point of g corresponds to the unstable 
periodic state of the oscillator. 

We will consider two major relaxation mechanisms of 
the oscillator: damping due to coupling to a thermal bath 
and dephasing due to oscillator frequency modulation by 
an external noise. Usually the most important damping 
mechanism is transitions between neighboring oscillator 
energy levels. They result from the coupling linear in 
the oscillator coordinate. Since the energy transfer is 
~ hujo, in the rotating frame the transitions look instan- 
taneous. We will assume that the correlation time of 
the noise that modulates the oscillator frequency is also 
short compared to l/l^wj, so that the noise is effectively 
i5-corrclated in slow time r. Then the quantum kinetic 
equation is Markovian in the rotating frame, 
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where Tp describes damping 

tp ^ r\duj\-'^ [{n+l){a'^ap-2apa^ + pa^a) 

+n{aa^ p — 2a) pa + paa})\ , (5) 

and fP^p describes dcphasing, 

pph^^pph|^^|-i ^~t^^ [a)a,p\] . (6) 

Here, T and PP^ are the damping and dephasing rates, 
a = (2A)~^/^((5 -I- iP) is the lowering operator, and n = 
[exp (hojo/kT) — is the oscillator Planck number. In 
what follows we use dimensionless parameters 



Q = |(5tj|/r, >?p^ = pp^/Ar. 
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We assume that >fi^ < 1. This means that the dephasing 
fluctuations intensity may be comparable to the intensity 
of quantum fluctuations associated with damping, which 
is oc Ar, see below, but that TP^ <C T. 

Metastable decay of the driven oscillator was studied 
earlier assuming that the damping-induced broaden- 
ing of quasienergy levels is small compared to the typi- 
cal interlevel distance. This condition necessarily breaks 



near a bifurcation point where local extrema of g{Q, P) 
come close to each other and the motion is slowed down. 
Therefore the analysis should be done differently. It is 
simplified in the Wigner representation of the density 
matrix, 



Pw{Q,P) = I d^e 



-'"'/'piQ + li^Q-lj], (8) 



where p{Qi,Q2) ~ {Qi\p\Q2) is the density matrix in 
the coordinate representation. Using Eqs. ©-Q one 
can formally write the equation for pw as a sum of terms 
proportional to different powers of A, 
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Here we introduced vectors K = {Kq,Kp) and V = 
(5q,9p). 

Vector K in Eq. Q determines the evolution of the 
density matrix in the absence of quantum and classical 
fluctuations, 

Kq = dpg - Q-^Q Kp = -dpg - Q-^P. (10) 

This evolution corresponds to classical motion 

Q = Kq, P = Kp. (11) 

The condition K = gives the values of Q,P at the 
stationary states of the oscillator in the rotating frame. 

The term L^^^ in Eq. Q describes classical and quan- 
tum fluctuations due to damping and dephasing. 
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These fluctuations lead to diffusion in (Q, P)-space, as 
seen from the structure of L^^\ 

The term L^^^ in Eq. describes quantum effects of 
motion of the isolated oscillator, 
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(Qdp - Pdg) V' 
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In contrast to L'--^\ the operator L'^' contains third 
derivatives. Generally the term X'^U^^pw is not small, 
because pw varies on distances ~ A. However, it becomes 
small close to bifurcation points, as shown below. 

From Eqs. (|10|l . (|ll|l . for given damping fi^^ the oscil- 
lator has two stable and one unstable stationary state in 
the rotating frame (periodic states of forced vibrations) 
in the range 

Pb\^1< P < Pb\^) and one stable state 



outside this range 
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At /S^' and /S^'' the stable states with large and small 
q2 _|_ p2^ respectively (large and small vibration ampli- 
tudes), merge with the saddle state (saddle- node bifur- 
cation). The values of Q,P at the bifurcation points 1, 
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In the absence of fluctuations dynamics of a classical 
system near a saddle-node bifurcation point is controlled 
by one slow variable |20l | . In our case it can be found by 
expanding Kq p in SQ = Q — Qb, SP = P — Pb, and 
the distance to the bifurcation point t] = P ~ Pb- The 
fimction Kp does not contain linear terms in SQ, SP. 
Then, from Eq. I|ll|l . P slowly varies in time for small 
6Q, 6P, rj. On the other hand 

Kq « -217-1 _ ^^^p) ^ ^ ^^^2Yb - 1). (16) 

Therefore the relaxation time of Q is il/2, it does not 
depend on the distance to the bifurcation point. As a 
consequence, Q follows P adiabatically, i.e., over time 
~ r2 it adjusts to the instantaneous value of P. 

The adiabatic approximation can be applied also to 
fluctuating systems. The approach is well known for clas- 
sical systems described by the Fokker-Planck equation 
[2H . We now extend it to the quantum problem. 

Formally we change in Eq. Q from Q and P to 5Q ~ 
5Q — qbSP and SP. For times r ^ the distribution 
pw has a narrow peak as a function of 6Q, whereas its 
dependence on SP is much more smooth. We seek pw 
near its maximum over SQ in the form 



27rAcr2)-i/2 pj.p fsQ^/2XaA pwiSP), (17) 
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where = i(l + a|) (n+ i) + I^^p^^Pb^^. The SQ- 
dependent factor in p^v is chosen so that in Eq. © the 
term dgqKqpy/ and the term oc X&^^^pw compensate 

each other. Note that corrections from A^L^^^pvK are of 
higher order in A for SQ^ < A. 

The function pw describes the distribution over SP. In 
the spirit of the adiabatic approximation, it can be cal- 
culated disregarding small fluctuations of Q, i.e., setting 
SQ = in Eq. 0. Formally, one obtains an equation for 
Pw by substituting Eq. (|17|) into the full kinetic equation 
^ and integrating over SQ. This gives 



Pw ~ dp [pwdpU + Wsdppw] 



where U and T) have the form 
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with b = -py'\2YB)-\l - 2n^YB + n^). In Eqs. (HHJ, 
(|19|l we kept only the lowest order terms in SP, P — Pb,X. 
In particular we dropped the term —X^QBdpPw/^ which 



comes from the operator L^^^ in Eq. One can show 
that, for typical \SP\ ~ 1??!"'^''^, this term leads to correc- 
tions ^ ?7, A to Pw- 

Eq. H18|l has a standard form of the equation for clas- 
sical diffusion in a potential U{SP), with diffusion co- 
efficient XVp. For rjb > the potential U has a min- 
imum and a maximum. They correspond to the stable 
and saddle states of the oscillator. The distribution pw 
has a diffusion-broadened peak at the stable state. Dif- 
fusion also leads to escape from the stable state, i.e., to 
metastable deca y. The decay rate W is given by the 
Kramers theory [23 , 



W = Ce"^^/^, Ra = 
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with prefactor C = 7:~^(br]/2)^/'^Pg^^^\SLu\ (in imscaled 
time t). 

The rate (|20|l displays activation dependence on the 
effective Planck constant A. The characteristic quantum 
activation energy Ra scales with the distance to the bi- 
furcation point -q — P — Pb as rf/"^. This scaling is inde- 
pendent of temperature. However, the factor Vb in Ra 
displays a characteristic T dependence. In the absence 
of dephasing we have Vb = 1/2Q. for n <Si 1, whereas 
Vb = kT/hw^VL for 1. In the latter case the expres- 
sion for W coincides with the result pll |. 

In the limit ^ 1 the activation energy H20(l for 
the small-amplitude state has the same form as in the 
range of P still close but further away from the bifurca- 
tion point, where the distance between quasienergy levels 
largely exceeds their width 0. We note that the rate 
of tunneling decay for this state is exponentially smaller; 
the tunneling exponent for constant quasienergy scales 
as ry^/* 10 , which is parametrically larger than rf/"^ for 
small T] [for comparison, for a particle in a cubic potential 
(|19|l the tunneling exponent in the strong-damping limit 
scales as 77 0]. 

For the large-amplitude state the quantum activation 
energy, Eq. (|2()|l . displays different scaling from that fur- 
ther away from the bifurcation point, where Ra oc /S^^^ 
for 3> 1 13 ■ For this state we therefore expect a 
scaling crossover to occur with varying p. 

The approach to decay of vibrational states can be 
extended to a parametrically modulated oscillator. The 
Hamiltonian of such an oscillator is 

H^{t) = ^ 1^2 [^2 ^Fco^iujFt)] + ^7^4 . (21) 

When the modulation frequency cjf is close to 2ljq, as 
a result of parametric resonance the oscillator may have 
two stable states of vibrations at frequency ujf/2 (period- 
two states) shifted in phase by tt [ijj. For <C the 
oscillator dynamics is characterized by the dimensionless 
frequency detuning p, effective Planck constant A, and 
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relaxation time C, 
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As before, A will be the small parameter of the theory. 

Parametric excitation requires that the modulation be 
sufficiently strong, ^ > 1. For such Q the bifurcation 
values of ji are 



(1.2) 
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= T(l-r')^/^ c>i. 
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If 7 > 0, as we assume, for /i < /i^'' the oscillator has 
one stable state; the vibration amplitude is zero. As /.t in- 
creases and reaches /x^"* this state becomes unstable and 
there emerge two stable period two states (a supercritical 
pitchfork bifurcation). They remain stable for larger /i. 

(2) 

In addition, when n reaches the zero-amplitude state 
also becomes stable (a subcritical pitchfork bifurcation). 
The case 7 < is described by replacing ^ — /i. 

The classical fluctuation-free dynamics for ^ close to 
/is is controlled by one slow variable The analysis 
analogous to that for the resonant case shows that, in the 
Wigner representation, fluctuations are described by one- 
dimensional diffusion in a potential, which in the present 
case is quartic in the slow variable. The probability W of 
switching between the period- two states for small ji — yL^g 

and the decay probability of the zero-amplitude state for 
f 21 

small 11 — iig have the form W — C exp(— i?^/A) with 
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The corresponding prefactors are 
'-2 



(2) 



I (1.2)\ 

(a*s = ) 

2C^^^ ^2^/'^^T-'^TC\^XB\\^l-^J.B\■ We note that dephasing 
does not affect the decay rate, to zeroth order in — pLg. 

From Eq. H24|l . at parametric resonance the quantum 
activation energy Ra scales with the distance to the bi- 
furcation point as 77^. In the limit 3> 1 the same ex- 
pression as Eq. (|24|l describes switching between period- 
two states still close but further away from the bifurca- 
tion point, where the distance between quasienergy levels 
largely exceeds their width. The exponent for tunneling 
decay in this case scales as 77 ^3 • 

It follows from the above results that, both for res- 
onant and parametric modulation, close to bifurcation 
points decay of metastable vibrational states occurs via 
quantum activation. It results from diffusion over a bar- 
rier. The quantum activation energy is smaller than the 
tunneling exponent. Near bifurcation points these quan- 
tities become parametrically different and scale as differ- 
ent powers of the distance to the bifurcation point. 

The exponent of the decay rate displays a characteris- 
tic dependence on temperature. In the absence of dephas- 
ing, for kT 3> ftwo we have standard thermal activation, 
Ra oc l/T. The low-temperature limit is described by 



the same expression with kT replaced by fta;o/2. Quan- 
tum activation imposes a limit on the sensitivity of bifur- 
cation ampliflers based on modulated Josephson oscilla- 
tors used for quantum measurements [^Q. 

In conclusion, we have studied decay of metastable 
states of forced vibrations of a quantum oscillator. Both 
energy dissipation from coupling to a bath and noise- 
induced dephasing were taken into account. We have 
found the exponent and the prefactor in the decay rate 
near bifurcation points. The quantum activation energy 
for resonantly excited period one states scales with the 
distance 77 to the bifurcation point as rf^'^, whereas for 
parametrically excited period two states it scales as 77^. 
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